The stationary solutions of G-P equations in double square well 
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^— ■») ' We present analytical stationary solutions for the Gross-Pitaevskii equation (GPE) of a Bose- 

I Einstein condensate (BECs) trapped in a double-well potential. These solutions are compared with 

those described by [Mahmud et al., PRA 66, 063607 (2002)]. In particular, we provide further 
evidence that symmetry preserving stationary solutions can be reduced to the eigenstates of the 
corresponding linear Schrodinger equation. Moreover, we have found that the symmetry breaking 
solutions can emerge not only from bifurcations, but also from isolated points in the chemical 
■ potential-nonlinear interaction diagram. We also have found that there are some moving nodes in 

the symmetry breaking solutions. 
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I. INTRODUCTION 



The Gross-Pitaevskii equation (GPE) describes many features of Bose- Einstein condensates (BECs) of dilute atomic 
gases in an external potential at zero temperature [H . The properties of the ground state of the GPE with external 
potential have been extensively studied; and many interesting phenomena have been reported. Of special interest 
is that the mean filed interaction profoundly modifies the density profiles and the stability of the ground state ^j. 
Recently, the properties of the nonCTound-state stationary solutions of the GPE have attracted more attention both 
theoretically and experimentally [Ijl, 01- For example, the dark solitons have been created in the atomic gases with 
^ , positive scattering length by phase engineering optical techniques i5|, Q . 
Q ' Based on Q, one can classify the stationary solutions of GPE as symmetry preserving and symmetry breaking 
O , solutions. Whereas the vortices and solitons observed in experiments [j, [E, are symmetry preserving solutions, 
which can be reduced to the eigenstates of the corresponding linear Schrodinger equations, the macroscopic quantum 
self-trapping state in the two states [§, [13, [HI, [l3| and non-Bloch state [l^ [iJI in periodical potential are symmetry 
J> i breaking solutions, and they can not be reduced to the eigenstates of the corresponding linear Schrodinger equations. 

I The system of GPE with double well is a good system to investigate the properties of the stationary solutions of 
\l GPE and also good for studying the special nonlinear dynamics, for example, the nonlinear self-trapping effect which 
was predicted theoretically in 1997 [Qj and realized experimentally [ll|, [3, [lB| last year. It has been shown that 
there are stationary solutions, which are either symmetry preserving or symmetry breaking, both numerically and 
analytically in this system 0, [^ . With the help of the double square well [2| , which allows one kind of analytical 
solution for GPE, Reinhardt and his collaborators have confirmed the numerical calculation for GPE with double-well 
trap [1| . This provides one possible way to investigate the properties of the stationary solution for GPE with double 
' well. It is interesting to know how the stationary solutions change as the nonlinear interaction increases or decreases, 
how the symmetry breaking stationary solutions emerge and how many kinds of the stationary solutions there are for 
I , ' the fixed nonlinear interaction. 

In this paper, we present new stationary analytical solutions for GP equation with double square well. Compared 
with the solutions of [3, ours can be reduced to the eigenstates of the corresponding linear Schrodinger equations. 
^ It is this feature which enables us to understand the above-mentioned interesting problems concerning the stationary 
solutions for GPE. It allows us to obtain the critical nonlinear interaction value, over which the symmetry breaking 
solution emerges. And from the relation of the profiles of the stationary solutions and chemical potentials with the 
nonlinear interaction, we can directly see the means of the symmetry and symmetry-breaking and different stationary 



$H ' solutions for special nonlinear interaction. 



II. MODEL AND THE SOLUTIONS 



Considering that the BECs of dilute atomic gases are confined by a very anisotropic harmonic potential (lo± 3> uj^, 
where uj±^x are the confined frequency in the y-z space and x direction) and that the BECs are loaded in a double 
square well in the weak confined direction (x direction) , then the dynamics of this system is governed by the 1-D G-P 
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equation 

^n^^ = + V (X) + 50 I* 1^ j * (X, t) (1) 

where go = ° is the 1-D reduced nonUnear interaction constant. The potential is of the form 

{oo |x| > a 
fe<|a:|<a (a = l/2,yo>0) 

Vo \x\ < b, 

The stationary solution can be written as ip {x) = r [x) exp (^—i^tj and after re-scaling the equations, we arrive at the 
equation for r (x) 
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fJ.r (x) = —-^r-^r (x) -+V {x)r (x) + -qr (x) (2) 

where 77 ~ Nogo ^"^f" and iVo is the total number of the atom. The energy and the potential are measured in unit of 

2^j^, and L is the length of the total space (here L — 2a). Due to the double well case, the stationary solution is 
just the real function, so we have assumed that our solution r {x) is real function. As in Ref. [2], the solution of ^ 
can be written in terms of the Jacobi elliptical function. 

Generally, we have two different solutions depending on the relation of the chemical potential and the barrier height. 
First for ^ > Vb 



ri{x)=A-S'N{Kx + 5,ni), (3) 

where 

"1 = 2^^, ti = K^ + Vo + —V- (4) 

This solution is also valid for Vq = 0, which corresponds to the region of 6 < < a. It is easy to check that 
when T] = 0, our solution is reduced to ri {x) — Asm {Kx + S) , where ni = 0, fj, = + Vq. This is nothing but the 
eigenstates of linear Shrodinger equations for /i > Vq. We have two different Jacobi functions for /i < Vb, corresponding 
to the region of |a;| < 6. To the case with one node in the barrier region, the solution is 

r2(a;) =B-SC(Q:e + 7,712), (5) 
n2 = l- f^ = Vo~Q'^ - —T]. (6) 



TABLE I: The modular transformation of Jacobi elliptic function 
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but to the case without node, 

r2(a;) =B-NC(Qx + 7,n2), (7) 



(8) 
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Same as the solution ([3]) in 77 = 0, our solutions are reduced to r2 (x) = i? sinh (Qx + 7) for ([5]) and r2 (x) = 
B cosh{Qx + for ([7]), where ^2 = 1, /x = Vq — Q^. It is interesting to note that those two solutions are precisely 
identical with the eigenstates of linear Shrodinger equations for one node or no node within the barrier. 

Please note here that we do not restrict the value of ni from to 1 as is usually used in Jacobi elliptic function. 
But this problem could be solved by the modular transformation table (see tableU where nii — j^^, m2 = ^.) p^ . 
The Jacobi elliptical function sc and nc are constructed from the Jacobi elliptical SN, CN and DN (see ftable HI)) or 

M)- 

III. SYMMETRY PRESERVING AND SYMMETRY BREAKING SOLUTIONS 

As mentioned in Introduction, we have two different kinds of stationary solutions depending on whether the sta- 
tionary solution has its linear counterpart @. The symmetry preserving solution has the linear counterpart, as it 
could be reduced to the eigenstates of the corresponding linear Schrodinger equations. But the symmetry breaking 
solution could not be reduced, therefore it does not have the linear counterpart. Such being the case, we are required 
to find the stationary solutions of Eq. ([T]). It is worthwhile to note that our method is just the usual one in linear 
case. 



TABLE IL Other Jacobi elliptic functions 



NS(u|m) 


_ 1 


SC{u 


m) 


_ SN(ii 


m) 




m) 


SN(u|m) 


— SN(u|m) 


— CN{u 


m) 


- DN(«|m) 


NC(?i|m) 


_ 1 


CS{u 


m) 


_ CN(u 


m) 


CD(u 


m) 


_ CN(u|m) 


— CN(u|m) 


— SN(ii 


m) 


~ DN(u|m) 


ND(«|m) 


_ 1 


DS{u 


m) 


_ DN{u 


m) 


DC(u 


m) 


_ UIN(u m) 


~ DN(u|m) 


— SN(-u 


m) 


— CN(u|m) 



With the help of ([3]), ([5|) and ([7]), solutions in the three regions can be written in the form 

fi{x) ^ Ai-S'N{Ki{x + a),ni), -a < x < -b 
f2{x) = S-NC(g(x + 7),m), \x\<b 
fsix) ^ A2-SN{K2{x-a),n2), b<x<a 

for the case without node inside the barrier. For the case with one node inside the barrier, we have 

fi{x) = Ai ■ SlSi {Ki{x + a),ni) , -a < x < ~b 
f2{x) = B-SGiQ{x + -,),m), \x\ < b 

fsix) = -A2 • SN(i^2(a; - a),^^2) . b<x<a 



(9) 



(10) 



We have considered the fact that the solutions vanish on and outside the potential |a;| > a. To fix the parameters 
Ai, A2, Ki, K2, Q, B and 7, we need the continuity 
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(11) 



and normalization conditions 



|/i (a;) I 



\fUx)\dx+ / \fi{x)\dx = l, 

b Jb 



(12) 



The definition of the chemical 



and we need that the chemical potential is the same in different regions {fii = fJ-2 — fJ-s) 
potential ^ can be found in Eqs.(j4]), Eqs.® and Eqs®. 

In Fig.l we present the first four symmetric and antisymmetric solutions for —100 < 77 < 100. It is clear that 
all of these solutions can be reduced to eigenstates of the corresponding linear Schrodinger equations (77 = 0). All 
of these solutions are symmetry preserving solutions. A barrier height of Vq = 1000, barrier width 2b — 0.1, well 
width 2a — I are used all through this paper. The shaded part represents the barrier region in all the figures. It is 
hence easy to see that the effect of the nonlinear interaction on the profile of the wave function of the high en ergy 
state is smaller than that in the case of the low energy state. This situation is the same as in the lattice case jl4| . 



FIG. 1: The symmetry preserving solutions of G-P equation 





FIG. 2; (Color online) The symmetry breaking solutions of G-P equation for positive nonlinear interaction 



In Fig. 4, the chemical potential have been plotted as the function of the nonlinear interaction. From Fig. 4, one can 
read that the chemical potential difference between the ground state and the first excited state increases with the 
increase of the strength of the nonlinear interaction. This seems to suggest that the positive nonlinear interaction 
enhances tunneling effect but the negative nonlinear interaction eliminates it. When the nonlinear interaction is very 
large and negative f r]^ —54.5), the symmetrical properties of the "ground" state change from the symmetry to the 
anti-symmetry 0, yj. In fact, when the nonlinear interaction is negative and large, the symmetry preserving solution 
is not the ground state. But this will be discussed in the later part of the paper. 

In Fig. 2 we present the first four symmetry breaking solutions for positive nonlinear interaction up to 100. When 
?/ = 100, these four solutions are exactly the one presented in In [H, they also show the profile of the first kind of 
solution in Fig. 2 for 77 = 15, 30, 50, 100. But here our results show that this solution emerges from the bifurcations 
of the chemical potential of the first excited state at ry^ ~ 0.32 (See Fig. 4). Some more detailed behavior can be read 
from Fig. 2. This solution has been predicted in [l3|. Its chemical potential is larger than that of the first excited 
state. It is easy to see that the node in this solution is moving as the nonlinear interaction increases. From 0.32 
to 5.563, the node is displacing from within the barrier to the right well. The fourth solution emerges from another 
bifurcation (77^ ^ 1-48), at which the chemical potential is the same as the fourth symmetry preserving solution 
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FIG. 3; (Color online) Same as Fig. 2 for negative nonlinear interaction 



(the third excited state). For this solution, there is one fixed node and two moving nodes. Due to the increase of 
the nonlinear interaction, one of the moving nodes moves from the barrier region into the well. Finally, there are 
two moving nodes in the right well and one fixed nodes in the left well. The second and third solutions in Fig. 2 
belong to a new kind of symmetry breaking solution. They occur from the same isolated point (77^ ~ 51.284) in the 
diagram of the chemical potential and the nonlinear interaction. Here an isolated point is one point whose chemical 
potential (or eigenvalue) is not the same as the neighboring points in the diagram of the chemical potential and 
nonlinear interaction, whether it is the symmetry preserving or symmetry breaking solution. Therefore their profiles 
are modified completely differently with the increase of the nonlinear interaction. The second one has one fixed node 
and one moving node while the third one has two moving nodes. 

The self-trapping effect as predicted in 1997 [l3| was experimentally realized in [HI, ESi [3 ■ This dynamical effect 
can be understood very well based on the two mode approximation [ifll . In our case, due to the zero phase difference 
between the two well wave functions, the condition for the self-trapping can be written as A = 2(1 -t- ^/l — K^)/ A^. It 
is easy to check that all of the symmetry breaking solutions satisfy this condition. We show this calculation in Table 
mil where Ac is the critical value for the self-trapping [lo| and A is our solutions. 



TABLE III: Different Numbers for the First Symmetric Breaking Solutions 
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5.0627 


0.4945 


0.9955 


10.0627 


0.4940 


0.9960 


20.0627 


0.4932 


0.9985 



In Fig. 3, the symmetry breaking solution for negative nonlinear interaction has been presented up to —100. When 
77 = —100, they are the same as the Fig. 4 in Q but with a different order. Same as the symmetry breaking solutions 
in positive nonlinear interaction, the first and fourth solutions emerge from the symmetry preserving solution at 
77c ^ —0.32 and at 77^ ~ —1.46 respectively. Please note that the first kind of symmetry breaking solution emerges 
from the ground state and its chemical potential is less than that of the ground state. Therefore this state is the 
ground state of the system (From Fig. 4). This is what is called the quantum phase transition in this system [2,^. 
This solution does not have any node and its profile show that the particle would stay in one of the wells for negative 
nonlinear interaction. The fourth solution includes one fixed node and a moving one. Now the node moves from the 
well into the barrier region with the decrease of the nonlinear interaction. Around 77 ~ —33.239 the node is in the 
barrier region. 

Again, the second and the third solution occur from the isolated point at 77^ ~ —43.08 in the diagram of the chemical 
potential with the nonlinear interaction. Now their nodes become fixed. Again the density profile has completely 
different behavior with the decrease of the nonlinear interaction. The second solution will collect the particles to one 
well but the third one keeps them in two wells. There is still another crossover between the third and the fourth 
solution near —100. That is why my order of the solution is different from at 77 = —100. 



6 




FIG. 4: (Color online) The chemical potential as the function of nonlinear interaction 



In Fig. 4. we plot the chemical potential as the function of nonlinear interaction both for symmetry preserving 
and breaking solutions. Due to the high barrier, the quantum tunneling effect between the two wells is small and 
the chemical potential difference between the ground state and the first excited state is small too. From Fig. 4, one 
can directly see that the symmetry preserving solution can be reduced to the eigenstate of the linear Schrodinger 
equations (77 = 0), while the symmetry breaking one can not. To see this more clearly, we amplify the figure around 
r] = (Fig(a) is for the ground state and the first excited state and Fig(b) is for the second excited state and the third 
excited state). The embedded figure in Fig. 4 (left) shows the critical point for the quantum phase transition [1,0]. 
We denote the symmetry preserving solutions as "L" , and the symmetry breaking solutions and positive nonlinear 
interaction as "NP" , and the symmetry breaking solutions and negative nonlinear interaction as "N" . It is clearly 
shown that the symmetry breaking solution can emerge from two kind of points. One is the bifurcation and the other 
one is the isolated point in this diagram. We arrange our solution in Fig.l, Fig. 2 and Fig. 3 in the chemical potential 
order. It is easy to see that the nature of the mean-field ground-state and the structure of the energy spectrum of the 
nonlinear system depend on the value of the nonlinear interaction. 



IV. CONCLUSION 



We have presented one new analytical solution which could be reduced to eigenstates of the corresponding linear 
Schrodinger equations. And we have shown the clear and direct evidence to the relation of the symmetry preserving 
stationary solution and the eigenstates of the linear Schrodinger equations. Based on the origin of the symmetry 
breaking state, one can find two kinds of different symmetry breaking solutions, neither of which can be reduced to 
the eigenstates of the corresponding linear Schrodinger equations and have one critical nonlinear interaction, over 
which the symmetry breaking solution emerge. One of the interesting things is that the first kind of symmetry 
breaking solution can exist under the following conditions: 1. for positive nonlinear interaction, it emerges from 
the first excited state and the third excited state; 2. for negative nonlinear interaction, it emerges from the ground 
state and the second excited state. We understand this based on the superposition of the eigenfunctions of the linear 
Schrodinger equations. And high accuracy of the critical value of the nonlinear interaction has been found for the 
first kind of symmetry breaking solution which emerges from the symmetry preserving one. This explains why we 
can not find the symmetry breaking solution from the ground state for the positive nonlinear interaction [l7| . As 
this method is just valid for the small nonlinear interaction, we can not yet understand the second kind of symmetry 
breaking solution. But how the profile of the stationary solution change with the nonlinear interaction has also been 
presented. Our calculation shows that the symmetry breaking solution satisfies the self-trapping condition. This will 
help to understand this effect from stationary solutions. 

Considering the situations realized in the experiment, our model can be regarded as too simple. But this model 
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can be solved analytically and is good enough to provide qualitative description. Usually, one can understand the 
realistic double well by the quartic function ax'^ + bx^ + c. The chemical potential or the eigenvalue of GPE may show 
a little difference from the results of our simple model, but they are consistent quantitatively 

WeiDong Li thanks the stimulative discussion with A. Smerzi and L. Pitaevskii. This work is supported by the 
NSF of China (No. 10444002), SRF for ROCS^EM of China and SRF for ROCS of the Shanxi Province. 

Note added. Recently, we noted two papers [isl. [l9l| which approach the same problem with similar techniques. In 
the present work we provide a more detailed analysis of the density profiles of stationary solutions as a function of 
the nonlinear interaction and we find a symmetry-breaking solution. 
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